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Abstract 

Island size and capture zone distributions (ISD, CZD) are studied numerically in submonolayer 
growth with various critical island sizes and shapes. CZD scaled by the variance show excellent 
agreement with the Wigner surmise, confirming the Pimpinelli-Einstein approach for large CZs / 
large island dynamics. The ISD decay as exp (—s'^), with 7 = 4, ~ 2.4 and 2 for point, fractal, and 
square islands, respectively. A scaling approach explains the values of 7 from the Gaussian decay 
of CZD and the efficiency of islands to capture diffusing adatoms. 
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Since island nucleation in the submonolayer regime determines several features of the 
subsequent growth of a thin film iH^, it motivated many theoretical works to explain the 
island size distributions (ISD) and the capture zone distributions (CZD). Theoretical distri- 
butions are frequently compared to experimental data from various materials and processes 
or with simulation data from atomistic models, providing estimates of parameters such as 
critical nucleus sizes and binding energies. Recent works on submonolayer growth of organic 
molecules js-T] and colloidal epitaxy js] increase the interest in the subject for the possibility 
of extending the knowlegde on atomic epitaxy. 

However, the exact forms of the CZD and ISD still remain unknown and, consequently, 
the basic mechanisms governing their kinetics are unclear. Most theoretical approaches 
assume irreversible aggregation of adatoms to islands whose size exceeds a critical value i. 
Important examples are that of Amar and Family (AF), who proposed an empirical formula 
for the ISD jol, the Mulheran and Blackman (MB) approach to relate ISD to distributions 
of areas of Voronoi polygons jlO|, and the recent proposal of Pimpinelli and Einstein (PE) 
ll| that CZD are described by the Wigner surmise (WS) from random matrix theory 12 |. 



The AF formula is widely used to fit the peaks of experimental ISD and to estimate critical 



island sizes (see e. g. Ref. [13|), but deviations from point island model distributions are 

nn 

significant [l^ . The solution of equations from the MB approach provides ISD very close 



to point island model simulations, but it is not so popular as the AF formula for fitting 
experimental data. The WS works better than the MB curves for fitting simulated CZD o " 
point and circular islands Uj , but recent numerical work by Amar and Evans groups 
showed deviations in the peaks and in the left tails. Alternatively, a Gamma distribution 
was used to fit CZD |17 |. 



Here we will analyze very accurate simulation data of ISD and CZD of growth models 
of point and extended (fractal and square) islands with irreversible aggregation {i = 1 and 
2). The CZD are fitted by the WS after rescaling by the variance, which reduces deviations 
in the peaks and left tails and hightlights the Gaussian right tail predicted by PE for all 
island shapes. The ISD shows deviations from the AF formula, but have universal right tail 
decays dependent on the island shape and related to the Gaussian CZD. This confirms the 
PE approach for describing the dynamics of large islands and large capture zones (CZs). 

We performed simulations in square lattices of size 2048 x 2048, and confirm negligible 
finite-size effects with runs in size 1024 x 1024. Diffusion-to-deposition ratios R = D/F 
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between 10^ and 10^° were analyzed, with critical island sizes i = 1 and 2 = 2. For i = 2, 
detachment probability e = 10^^ was used with point islands to facilitate island growth and 
reduce fluctuations, but for extended islands we used e = 1. Averages are taken from 10^ 
lattice configurations for each i, R and coverage 6. In point island models, all atoms of an 
island a single lattice site. The model version 

was simulated, with an adatom aggregating to an island as it is deposited or hops onto it 
(however, simulations of the neighbor-attachment version 1^ leads to similar conclusions). 
In fractal island models, each atom permanently aggregates at the site where it collides 



with an island 20| , ^nerating branched clusters that resemble those of diffusion- limited 
aggregation (DLA) [21,]. In the square island model, instantaneous relaxation to compact 
(square) shape is perfomed after an adatom collides with an island {2^. A site belongs to 
the CZ of an island if its distance to the island border is smaller than its distance to any 
other island border. The CZ area x is the number of sites in the CZ, and the island size s 
is its number of atoms. 

It is expected that the probability density of CZ area x follows the scaling form 

and an equivalent relation is assumed for the density of islands of size s, Q{s). The alternative 
scaling of CZD with the variance ax = {{x — {x)) ) as 

P{x) = -g{u) , u^^, (2) 



was used in a recent study of roughness distributions l23\ to reduce the effect of scaling 
corrections in x and ax (here, corrections are expected to arise from small CZs, whose 
dynamics continuous theories fail to describe). For comparison, the WS is 

Pl3{x) = apx^ exp , (3) 

with the parameter /3 = | (i + 1) predicted by the PE theory, where d is the substrate 
dimension {d = 2 m this work). The AF formula for ISD is 

fi{u) = Ciu' exp (-miM^/''' ) , (4) 

where u = s/ (s) , Ci and aj are normalization constants. 
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In Fig. la we show scaled CZD for point islands with i = 1. An excellent collapse of 
the data for different R is observed with scaling by the variance, as shown in the inset of 
Fig. la. This is an important starting point for comparison with (/^-independent) theo- 
retical approaches. This contrasts to the scaling by the average, which shows significant 
/^-dependence of the peak heights [15]]. Fig. la also shows four theoretical curves for com- 
parison with the simulation data. The WS with (3 = 2, predicted by PE theory, is a good 



fit for the peaks under this resca 
previously observed in Refs. 15 



ing (inset of Fig. la)), but shows deviations in the tails, as 



16| . The Gamma distribution that best fits the CZD peaks 



has remarkable differences in both tails. The WS with /3 = 3 was suggested in Refs. |16L 124 1 
and turns out to be a good fit for the peak and left tail of the simulation data, but also 
differs from the right one. Finally, the generalized Gamma (GG) fit g{x) ~ x'^ exp {—bs"') 



with /3 = 4 and n = 1.5, proposed in Ref. |l6j, has small deviation only in the right tail. 

Fig. lb shows the same simulation data of Fig. la with abscissa u"^. The fit for u"^ > 15 
highlights the Gaussian right tail. This is the universal decay of CZD predicted by PE 
theory. 

Fig. Ic shows scaled CZD for point islands with i = 2, compared to the WS with /3 = 3, 
as predicted by PE theory. Good data collapse is also obtained with this rescaling (inset of 
Fig. Ic). There is discrepancy from the PE curve in the left tail, but the agreement in the 
right tail is very good, which confirms the Gaussian decay. 

The comparison of the skewnesses ((x — {x))^) / cr^ of the simulated CZD and the WS gives 
additional support to the PE approach: simulation gives S = 0.51 ±0.01 and S = 0.41 ±0.01 
for i = 1 and i = 2, respectively, while the values of the WS are 0.486 for /3 = 2 and 0.406 
for P = 3. It is also observed that CZD for point islands have a negligible dependence on 
the substrate coverage 6 when scaled by the variance. 

Scaled CZD of fractal islands with i = 1 are shown in Fig. 2a, with the predicted WS 
curve (/3 = 2). The good fit in the tails confirms the Gaussian decay, and scaling by the 
variance also leads to a good data collapse for different R and agreement with the WS peaks. 
Results for square islands with i = 2 are shown in Fig. 2b, again in good agreement with 
the respective WS curve (/3 = 3). Similar good fits by WS curves are obtained for fractal 
islands with i = 2 and square islands with i = 1. 

For fractal islands, the good data collapse shown in Fig. 2a is obtained for coverages 
9 < 0.20. For square islands, sub-island overlap enhances coalescence effects, thus the 
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behavior shown in Fig. 2b is observed only for 9 ^ 0.10. Also note that, for small R 
(typically R ^ 10^), the right tails for extended islands deviate to simple exponential decays, 
and data for i? = 10^ show crossover behavior. 

These results are striking evidence that the PE theory actually captures the main ingre- 
dients of the dynamics of large CZs. Indeed, the phenomenological model introduced by 
PE fl]]| extracts the Wigner distribution from a Langevin equation, which is a continuous 
approach expected to apply for large x. The deviations in scaling by the average (Eq. [1]) and 
the better data collapse in scaling by the variance (Eq. [2]) confirm that scaling corrections 
are reduced by the latter. Those corrections are related to small CZs, whose dynamics is 
not accurately described by the continuous PE approach. For instance, this explains the 
larger deviations f or p oint islands with z = 1, particularly in the left tails (as previously 



observed in Refs. 



15 



16|): since adatoms have to move for long times before reaching an 



existing island, the probability of two diffusing adatoms to meet (nucleating a new stable 
island) is large, leading to division of large CZs into smaller ones. For i = 2, the instability 
of two- atom islands reduce this effect. 

The scaling by the variance makes the peaks of CZDs indistinguishable for different i?s 
in a linear plot, in contrast to the log-linear plots that highlight the differences in the tails. 
Thus, it will probably be not useful to distinguish the critical nucleus in experimental data, 
where ?uctuations in the tails are large. However, it is certainly important for comparisons 
in theoretical works. 

Now we present results for ISD. 

Fig. 3a shows results for point islands. The variable (s/(s))^ highlights the right tail 
decay as exp (— s^) for large R. The inset of Fig. 3a shows the same data scaled by the 
variance as [with z = {s — {s))/crs]. This rescaling does not provide collapse of ISD data for 
different R and 6 and introduces corrections in the right tail decay shown in the main plot. 
Fig. 3a also shows the AF formula, which does not represent these ISD because it gives 
Q{s) — 7- as s/ (s) — 7- 0, as noted before 

For fractal islands, scaling by the variance provides good data collapse for different R, as 
shown in Fig. 3b. The AF formula is a good fit for the ISD peaks (inset of Fig. 3b), which 
explains its wide use, but it does not represent the right tail (main plot). The decay of the 
ISD is only slightly faster than the Gaussian decay of the CZD. The good fits of the right 
tails by exp {—z'^-^) (main plot) or exp (— s^ '') (if scaled by the average) explain the choice 
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of the abscissa in Fig. 3b. 

The ISD for square islands are shown in Fig. 3c. Good data collapse for different R is 
obtained in the right tails with scaling by the variance, despite the large differences in the 
left tails and peak heights. The same Gaussian decay of the CZD is obtained. Again, the 
AF formula shows deviations from the simulated ISD. 

The ISD for extended islands (fractal or square) have significant coverage dependence. 
For very low coverage, point island decay [~ e(~* )] is found, and, for large coverage, the 
right tails tend to simple exponentials. The behavior shown in Figs. 3a-c is typical of a 
narrow (though widely studied) range 0.05 < 9 < 0.2. Analogous results are obtained for 
i = 2. 

Now we establish a connection between the CZD and the ISD which explains the observed 
decays of the latter. Since the decay of CZD is universal, as predicted by PE, it is natural 
to expect that the decay of ISD will be related to that Gaussian decay and to the particular 
island shape. Working with large islands and capture zones is certainly necessary for such 
continuous approach. 

For fractal islands, it is plausible that larger islands have larger CZ. This is shown for 
circular islands in Ref. [25]. Thus, consider that islands of size s have CZ of typical size x, 
so that P{x) ~ Q{s). The shape of that island resembles that of a DLA cluster for small 
coverage because most of the attached atoms came from the island neighborhood by diffusion, 
and not by deposition inside that island. Thus, the typical radius of the island is i? ~ s^^^^, 
with the fractal dimension Dp ^ 1.694 of DLA 26|. The area of the CZ is expected to scale 
as the area of a circle with that radius, thus x ~ i?^ ~ s^^^^ . Consequently, the Gaussian 
tail of the CZD implies a tail for the ISD as exp (— x^) ~ exp s^^^^^ exp (— s^''^^). This 
faster-than-Gaussian decay is in excellent agreement with the results in Fig. 3b. 

For compact (square or circular) islands, the fractal dimension is Dp = 2, thus the 
same arguments lead to a Gaussian decay of the ISD, again in excellent agreement with our 
numerical results (Fig. 3c). 

For point islands, consider a typical CZ of area x of an island with size s {s atoms in a 
single site), with x, s ^ 1. The radius of this CZ is i? ~ x^/^, and its border has length of 
order R. It is surrounded by other CZs, whose number grows proportional to the border 
length R (typically they are small CZs neighboring a large one).. Now consider islands whose 
CZ is of size x + Ax and size s + As. The increase in the size is due to atoms deposited in 
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the additional area Ax which were able to reach this point island, instead of being absorbed 
by the neighboring islands. We expect that a fraction ~ of those deposited atoms will 
actually reach the central point island. Thus, As ~ '^^^^ which gives ^ ~ -Jj2 in the 
continuum limit and, consequently, x ~ (this means that large CZ area x corresponds to 
not so large island size s ~ x^/^). The Gaussian tail of the CZD implies a tail for the ISD 
as exp (— x^) ~ exp (— s^). This is also in excellent agreement with our numerical results for 
i = 1 (Fig. 3a) and reasonable agreement for i = 2. 

In summary, we studied numerically the ISD and CZD in submonolayer growth of islands 
with different shapes (point, fractal, and square) and critical sizes i = 1 and i = 2. The 
Wigner surmise describes the CZD when scaled by the variance of CZ area, with deviations 
in the left tails for point islands with i = 1. This confirms the super universal Gaussian 
decay predicted by Pimpinelli and Einstein, in contrast with the fits proposed by other 
approaches, and improves recent numerical work which did not focus on the CZD decay 
and used scaling by the average. The ISD are also not represented by well known fitting 
formulas, but universal decays of the right tails are found for each island shape and are 
explained by connections to the Gaussian tails of the CZD. We believe that these results 
motivate experimental efforts to measure accurate ISD and CZD, which are possible with the 
advance in microscopy techniques It also motivates theoretical efforts to understand 
ISD features other than the right tails, the reversible growth processes of atomic islands, 
and growth of more complex nanostructures Q]. 
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FIG. 1: (Color online) (a) Scaled CZD for point islands with i = 1 and ^ = 0.10 [R = 10^ (crosses, 
red), R = 10^ (triangles, blue), and R = 10^ (squares, green)] compared with the WS with /? = 2 
(solid curve), the Gamma distribution (dashed curve), the WS with /3 = 3 (dotted curve), and the 
GG fit with /3 = 4 and n = 1.5 (dashed-dotted curve). The inset is a linear-linear plot of the same 
data and the WS with /3 = 2. (b) Same simulation data of (a) with abscissa . The solid straight 
line highlights the Gaussian tail, (c) Scaled CZD for point islands with i = 2 and the same values 
of 9 and R, compared with the WS with /3 = 3 (solid curve). 
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FIG. 2: (Color online) Scaled CZD for (a) fractal island model with i = 1 and (b) square island 
model with i = 2, with R = lO'^ (crosses, red), R = 10^ (triangles, blue), R = 10^ (squares, green) 
and R = 10^^ (stars, magenta). The solid curves are the WS with (a) /? = 2 and (b) /3 = 3. The 
insets are linear-linear plots of the same data. 
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FIG. 3: (Color online) Scaled ISD for (a) point islands, (b) fractal islands, and (c) square islands, 
with 1 = 1 and 9 = 0.10. Symbols correspond to i? = 10^ (crosses, red), R = 10^ (triangles, blue), 
R = 10^ (squares, green) and R = 10^*^ (stars, magenta). The insets show linear-linear plots of the 
same quantities. The solid curve in all plots is the AF formula with i = 1. The dashed straight 
lines in the main plots are guides for the eyes. 
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